Circular Orbits in Einstein-Gauss-Bonnet Gravity 



OO 
O 

o 



O 



> 
1^ 



Valeria M. Rosa* 

Departamento de Matemdtica, Universidade Federal de Vigosa, 36570-000 Vigosa, M.G., Brazil 

Patricio S. Letelier^ 

Departamento de Matemdtica Aplicada-IMECC, Universidade Estadual de Campinas, 13081-970 Campinas, S.P., Brazil 
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I. INTRODUCTION 



The physical laws for nonrelativistic phenomena are usually written in a three dimensional space and for relativistic 
phenomena in a four dimensional spacetime. In both cases these laws admit natural extensions to spaces with arbitrary 
— ' number of dimensions, these can be taken as the possible existence of an operative principle (or principles) which in 
additions to these laws lead us to determine the "true dimensions" of the spacetime, the arena where the physical 
phenomena occur. 

Kant observed that the three-dimensionality of the usual space could be related with Newton's inverse square law; " 
. , The reason for the three-dimensionality of the space is yet unknown. It is likely that the three-dimensionality of space 
^jrj" results from the law according to which the forces on the substances act into each other. The three-dimensionality 
seems to result from the fact that the substances of the existing work act into each other in such a way that the 
strength of the effect behaves like the inverse of the square of the distances" [l| ■ 

Since in a first approximation, the galaxies and planetary systems can be considered as formed by particles moving 
in circular orbits around an attraction center, a possible principle to determine the dimensionality of the space is 

^ ' the existence of stable circular orbits. The study of the stability of circular orbits moving in potentials solutions of 

the n-dimensional Laplace's equations can be found in [sj. These studies on stability culminate with the Ehrenfest 
work [J] where he shows that in the Newtonian-Keplerian problem there exists stable limited orbit if and only if the 
dimension of the space is three (the potential must be null at infinity) . 

■ In Einstein's gravity Tangherlini [5)] study the stability under radial perturbation of orbits associated to test par- 
00 \ tides moving around a n-dimensional Schwarzschild black hole. He found stability for n = A only, i.e., in complete 

■ concordance with the Newtonian case. The dimension of the constant time slice, n — 1, will be denoted by d. 
Generalization of the Newton theory that consists in the addition of a cosmological constant was studied by several 

authors, see for instance [s] (cosmological context) and (potential theory context). Clearly this generalization is 
' motivated by the success of Einstein's gravity with cosmological constant. An example of model that makes use of 
a cosmological constant is the A Cold Dark Matter model (ACDM) S] that explains in a simple way the observed 
acceleration of the Universe as well as the cosmic microwave background. 

The Hilbert Lagrangian (Ricci scalar times square root of the metric determinant) is a divergence in two dimensions. 
We have that the Hilbert action is proportional to the Eulcr characteristic of the two dimensional manifold (Gauss- 
Bonnet theorem Q ) . In four dimensions we have that the Bach-Lanczos Lagrangian [HI (a particular combination 
of terms quadratic in the Riemann-Christoffel curvature tensor) is a divergence, so in four dimensions this Lagrangian 
has no dynamical significance, like the Hilbert Lagrangian in two dimensions. The simplest generalization of the 
Hilbert Lagrangian for n-dimensional spaces (n > 4) can be obtained adding the Bach-Lanczos Lagrangian to the 
Hilbert one. This new theoryj)f gravitation, Einstein-Gauss-Bonnet (EGB) theory, has some remarkable mathematical 
properties, see for instance 12']. 

Another outstanding feature of the Einstein-Gauss-Bonnet theory is that the Bach-Lanczos action appears as the 
leading quadratic correction term of the Hilbert action in the expansion of the supersymmetric string theory (isj . To 
be more precise the Bach-Lanczos action appears exactly in the expansion of the Eg, x heterotic string model [l3| • 
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The aim of this paper is to study the stabihty under radial perturbations as well as vertical perturbations of circular 
orbits related to test bodies orbiting a spherically symmetric center of attraction in the above mentioned theories of 
gravity in an arbitrary number of dimensions. 

In Section II we study the stability of circular orbits in a rt-dimensional Newtonian context with cosmological 
constant. In Section III we consider the stability of circular orbits in the EGB theory. We divide the study of stability 
in two. First, in Sect. IV, we consider the case without cosmological constant. In the next section (Sect. V) we 
analyze the case where this constant is not null. We end summarizing our results in Sect. V. 



II. d-DIMENSIONAL NEWTONIAN CENTER OF ATTRACTION WITH COSMOLOGICAL 

CONSTANT 



The Laplace equation with cosmological constant is 

V2$ + A = 0, (1) 

where is the Laplacian operator in d-dimensions, $ and A are the Newtonian potential and the cosmological 
constant, respectively. For spherical symmetry Eq. ([1]) reduces to 

that has as solution for n > 4 , 

CIA, 

$ = r 

n-3r"-3 2(n-l) ' ^ ' 

where C is an integration constant, that for n = 4 (d = 3) it is proportional to the mass of the center of attraction 
(C = GM, with G the Newton's constant of gravitation). 

The Newton equations for an axially symmetric potential can be written as 

«-«.'^-||, (4) 
Jt^"'^' = 0. (5) 

where the overdots denote derivation with respect to t, {R, (f, z'^, z"*"^) are d-dimensional cylindrical coordinates, 
and $ = z^, Defining the effective potential <I>e// = <i> + /i^/(2i?^), with h = i?^<y9(=constant), we find, 

dR ' dz^ ■ ^ ' 

Now let us study the stability of the circular orbit R — Rq, z'' — Q. We shall assume that the potential has reflection 
symmetry on the planes z' = 0, i.e., $ — ^[R, (z^)^, (z''^^)^]. Expanding the effective potential around this circular 
orbit we find, 

* (R 0) I l g^'i>e//(fi0,0,...,0) .„ , li^ d^^,ff{Ro,0,...,0) , 

'Pe// - 'Pe//(-Ko,0, ■-,0) + ^ [R~Ro) +2/^ ^77)2 > +■■■• 

From Newton's equations we find the motion for the perturbations p = R — Rq and = — 0, 

P + K^p = 0, d+^!C ^0, (9) 

where 

,.2 _ ^'^e//(-Ro,0,-,0) _ g'$(i?o,0,-,0) , 3 9^(flo,0,...,0) 

dR^ ^ dR^ ^ R dR ' ^ ' 

diz'Y d{z^y ' ^ ' 
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Thus when > we have that the radial perturbation remains bounded, i.e., we have hnear stabihty. This 
perturbations are usually named as epicyclic perturbations (T5| . In a similar way when vf > we have that the 
"vertical" perturbation also remains bounded. In this case we say that we have linear stability under perturbations 
along the direction z\ For the spherically symmetric potential ^ we find 

= (5-n)Ci?^-"-4A/(n- 1), (12) 
1/2 = ,y^ = CRl,-" - A/(n - 1). (13) 

To derive the previous equations we used r = [R"^ +'l2'i=ii^^)'^]^ ■ Note that because of the spherical symmetry of the 
potential we have that all the frequencies for all the vertical perturbations are equal. 

The relations and (fTS]) tell us that when A = we have stability only for n ~ 4. For A > and n = 4 we have 
stability whenever Rq < 3C/(4A). These last relation puts an upper limit to the size of the stable galactic structures. 
For A > and n > 5 we have < and the circular orbits are not stable under epicyclic perturbations. Now for 
A < we have stable circular orbits for n — 4,5 and for n > 6 we have stability for Rq limited by the lower bound 
Rq = [{n — 5)(n — l)C/(— 4A)]-'^/("'~^)(< Rq). In other words for certain values of n the presence of a cosmological 
constant can stabilize the circular orbits, sometimes with limitations in the size of the orbit. 



III. EINSTEIN-GAUSS-BONNET THEORY OF GRAVITATION 



The Einstcin-Gauss-Bonnct action is. 



S = / (Tx^ 



■\{R-2A + aLGB) 



'matter' (14) 



where R is the n-dimensional Ricci scalar (n > 4), k„ is defined as k„ = y/SnGn, where G„ is the n-dimensional 
gravitational constant. We use units such that the speed of light is one (c = 1). The Gauss-Bonnet term Lqb is given 
by 

Lgb = i?' - ^R^.R^" + R^.psR^""', (15) 

and a is a coupling constant. -Smatter P^^^ action that describes the matter. We use the conventions, 

^bcd — '^c^bd + ■ ■ ■ ) Rfj.!^ = ^rrv ^P^^ R = Rli- The action (fH|) is the simplest one for n > 4 built with topological 
terms. For the general case see p^ . 
From ([HI) we find 

Gp,u - aH^i, - Kg^,^ = i^nTf_cu, (16) 

where 

Hfj,i, = 2 [RRf_ii, — 2R^aRZ — 2R°''^ R^iaiyp + R"^'^ RuafSf] — -^^giiuLcB- (17) 

For vacuum (T^i/ = 0) spherically symmetric solutions to the EGB equations there is not lost of generality in 
choosing the metric, 

ds^ - f{r)de ~ f{r)-^dr^ - r^d^l^^, (18) 

where d^^_2 — '^^'i + ^j=3 (nr=r2^ Qi)d9^^ is the metric of a (n — 2)-sphere. The case in which this sphere is changed 
by a space of negative or null curvature has also been considered [l8[ . 

From the fact that the the Einstein tensor, G'"', as well as the Gauss Bonnet tensor, H^^ , are divergence free [lB| 
we have that equation (fTB)) tells us that T^^ is also divergence free. From this last fact we conclude that test particles 
move along geodesies [17i] . We shall analyze the stability of the circular geodesies 7 whose parametric equation is, 

t = t{s), r = ro, J = 2,...,n-1, 0^ ^ cj) E [Q,2n]. (19) 

The evolution of a small perturbation, of 7 is given by the geodesic deviation equation, 

^ + RLp^''ru' = Q^ (20) 
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where D/ds = u^\7 ^, and = dx^ /ds. Equation (|20p can be cast in a simpler, but not manifestly covariant way, as 



where the comma denotes ordinary differentiation. For the curve (|19p the system (PT|) reduces to, 

|i + n^e = 0, 

+ = 0, 

|n + 1^^1=0, 



(21) 

(22) 

(23) 
(24) 

(25) 



where 



f{ro){rof"{ro) + 3f (rp)) - 2rof (rp)^ 
ro(2/(rp)-rp/'(rp)) 
f'iro) 



ro(2/(rp)-rp/'(rp)) 



(26) 



We shall do a study of cases depending on the values of a and A. The analysis starts with the stability of circular 
orbits for the pure Einsteinian case (no Gauss-Bonnet term). Then we study the stability of these orbits when a 
coupling term and/or the cosmological constant are considered. 

For A = and a = 0, the model of spacetime reduces to the n-dimensional Schwarzschild solution. The equation 
for the metric function is 



r/'-(n-3)(l-/) = 0. 



(27) 



which has the general solution /(r) = 1 — C/r 



C is an integration constant that will taken as positive in order 



to have the Schwarzschild solution when n = 4 (C = 2GM). In this case the constants i/^ and can be cast as, 

C(n-3)ro"-3 



n — I/O TL—3 



C{n-1)) 



C{n - 3)((5 - n)r^-3 _ (^(^ _ ^^^^ 



r^^(2r^--C(r^-l)) 

Therefore, we have i/^ > if and only if r'^^^ > C(n-l), this imply that > if and only if (5-n)ro"^-C(n-l) > 
0, and it only occurs when n = 4 and rp > 6GM. Hence, we have stable circular orbits only in a 4-dimensional 
spacetime as it was expected. 



IV. EINSTEIN-GAUSS-BONNET THEORY WITH A = 0. 

Taking A = in with T^,^ = 0, the equation of metric function / is 



r/'-(n-3)(l-/) 
which has the general solution 



1 



a{n ~ 4) 



(2r/'-(n-5)(l-/)) 



0, 



/±(0 = i + 



2a 



l± \ I 



— 1 



(28) 



(29) 



where a = (n — 3)(n — 4)a and C is an integration constant. The sign of C is chosen analyzing the behavior of / for 
|a| 0. Let /+ be the function / where the sign of square root is positive and let /_ be the other case. When a is 
small, the asymptotic behavior of /_ is given by 



(30) 
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FIG. 1: GR branch, A = case. Epicyclic and vertical frequencies for several values of n and C = 1, 1 < ro < 3. For the first 
two figures a — 20 and for the third a = 80. In the first graphic n take the values n = 9 (top curve), 8,7,6 and 5 (bottom 
curve). In the second, counting along the line ro = 2, n takes the values 5 (top curve), 6, 7, 8 and 9 (bottom curve). And for 
the third graphic, counting along the same line, n — 6 (top line), 5, 7, 8, 9, 10 (bottom line). We see that the region of stability 
shrinks and disappears for large n. 



So, considering the Gauss-Bonnet term as a perturbation of the n-dimensional Schwarzschild solution, C must be a 
positive constant. 

When a is small, the asymptotic behavior of /+ is given by 

For a > 0, it corresponds to Schwarzschild-anti-de Sitter spacetime with negative gravitational mass, with the standard 
energy definition in background anti-de Sitter p^ . 

Then, there are two families of solutions corresponding to the sign in front of the square root that appears in (I29p . 
Following the definition given in [20*1 , the family with minus (plus) sign will be called general relativity (GR) branch 
(non-GR branch) solution. The spacetime structure of these solutions of the EGB theory is studies in [21i] . 



A. The GR branch solution. 



By chosing the solution / = /- of ((28|) from ([26]) we find. 



rj(l- V^) + Cro5(5-n) 



K 



,2 



-(r3C(n - 1) - 2r^VP) 
A1+A2+ A3 



-2r'^P^a[^2r^^ + rlC{n - 1)) ' 
where /3 = 1 + AaCrl~" and 



(32) 



A, = 8r2"(l - V^) + r^^+^C{n' - 1)(1 - 

A2 ^ 2rJ,'+^[32 - y^in^ - 8n + 39)] + 2^0^ a{n - l)[n + 7 - ./f3{ji + 3)], 

A3 = ArlC^a'^[32 - 2^aC{n - 5)(n - 9)] - 8r^""Ca(n - 5)(n - 1). (33) 

The behavior of the functions and are not simple, but it is possible to find n, a, C and rg such that and 
are both positive, i.e., there are stable circular orbits. For example, we have stable circular orbits when n — b, 
5 = 2, C = 1 and 0.2 < ro < 2.0; n = 6, a = 20, C = 1 and 0.2 < ro < 2.0; n = 7, <5 = 20, C = 1 and 0.7 < ro < 1.9; 
n = 8, <5 = 20, C = 1 and 1.07 < ro < 1.47; and n = 9, a = 80, C = 1 and 1.33 < ro < 1.46. In all these examples we 
have /_(ro) > in the specified region. Hence, the Gauss-Bonnet term can stabilize circular orbits in some higher 
dimensions spacetimes. 
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FIG. 2: Non-GR branch, A = case. Graphics of the functions v'^ and for a = 20, C = 1 and 1 < ro < 7 and n = 10 (top hne), 
9, 8, 7, 6, 5 (bottom hne). The have a large region of stabiUty. 

In Fig. [T]we show the epicychc and vertical frequencies for several values of n and C = l, l<ro<3. For the first 
two figures a = 20 and for the third a — 80. In the first graphic n take the values n — % (top curve), 8, 7, 6 and 5 
(bottom curve). In the second graphic, counting along the line = 2, n takes the values 5 (top curve), 6,7,8 and 
9 (bottom curve). There is no stable circular orbit when n = 9. And for the third graphic, counting along the same 
line, n = 6 (top line), 5,7,8,9, 10 (bottom line). In this case here is no stable circular orbit when n = 10, but for 
n = 9 we have stability. 

The vertical frequency, v'^ , is always positive no matter the value of n but the range of tq where is positive 
decreases when the dimension n increases. 

For n > 10 and a > 0, and also for n > 5 and ci < 0, we calculate the functions and into a grid of points with 
0.1 < r-o < 2 (step 0.1), 5 < n < 20 and 0.01 < a < 100 (step 0.01 into [0.01, 0.1], step 0.1 into [0.1, 1], step 1 into 
[1, 100]) and C ~ 1. In the case a < we work with —100 < a < —0.01, with the same steps for the corresponding 
intervals. In all these numerical examples we obtain ly^K^ < and, therefore, no circular orbit is stable. 

B. The non-GR branch solution. 

For completeness we study the case when / = /+ is the solution of In this case the constants i'^ and are 

given as in ([5^ and ([55]) replacing —y/P by y/j3. 

Also in this case, the stability depends on the sign of a. When a < 0, in the expression for — N^/Di,, we have 
a numerator, = Tq [(1 -I- ^fP)r1^ — [n— 5)Caro], non negative and a denominator, = a{rQC{n — 1) + 2rg v^), 
always negative. So, v'^ is not positive, hence there are no stable circular orbits. 

The parameters n{> 5) and C are positive. When a > the denominators of and are always positive 
(£)„ = 2(3^/Pr'^D^). The same happens with the functions /+, /?, and with the coefficients of tq in Ai and A2 (Note 
that Ai is the function Ai after replacing — by \fp, i = 1, 2, 3). These positive coefficients assure that, for ro 
sufficiently large, the numerator of , = [{l + ^/]3)rQ — {n — 5)CarQ], and the numerator of , A^k = [^1+^2 + ^3], 
are positive. Hence, we conclude that it is possible to find stable circular orbits for n > 5 and a > 0. 

Fig. [5] shows the functions 1/^ and , for C = 1, a = 20, 1 < r < 7 and n = 10 (top line), 9, 8, 7, 6, 5 (bottom 
line). We have stability for tq > r„, where r„ depends on the dimension n. 

V. EINSTEIN-GAUSS-BONNET THEORY WITH A / 0. 

Taking a = and A 7^ in (jl6p with = 0, the metric function obeys the equation, 

2Ar2 

r/'-(n-3)(l-/) + - = 0, (34) 

n — 2 
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which has the general solution, 



C 2Ar^ 

/W = 1 - - („,i)(^„2) ' 



the integration constant C will be considered positive as before. The frequencies ([^5]) in this case are = Ni/Di 
and = N2/D2, where 

Ni = Cr2(n - l)(n - 2)(n - 3)^ - 4Ar^'+i; Z^i = r;;(n - l){n - 2)(2rJ-3 _ _ ^^^^^ 
iV2 = Cr2(n - 2)(n - 3)2((5 - n)r[f-^ - C(n - 1)) - 2Arl^+\8ri;-^ - Cin^ - 1)), 

D2 = r'o^\2r'^-^ - C{n - l))(n - l){n - 2). (36) 



When A < 0, the numerator Ni of i/^ is positive for all tq, and its denominator Di is positive if and only if 



„n— 3 



> C{n — l)/2. Note that D2 change of sign with Di. The function / is positive when rg~ > C, and this 
condition is satisfied when Di > 0. So, it is necessary to find the values of tq for which the numerator iV2 of is 
positive. For better visualization, we write A = —A in N2 



N2 = 16Arg"-2 



2C A(n' - l)?-o - C{n - 2){n - 3)(r 



C2(n-2)(n-3)(7i-l)rg. 



Analyzing the variation of the signs of the polynomial N2, we can see that there is at least one positive root for 
A^2- Adding the fact that for large ro, ^^"2 is positive, we can conclude that it is possible to find a value f such that 
> Q and k? > Q for > r. Therefore, when A < and for n > 4 there exist stable circular orbits. 



When A > 0, in order to have / > it is necessary that, at least, 

2r"~3 „ 0(71-^1) > 0. 



(37) 



Under this condition the denominators Di and D2 are positive. Now it is necessary that the numerators A^i and N2 
be also positive. The polynomial A'^i is positive when 



Cr^(n- l)(n- 2)(n- 3)^ 



4A: 



> 0. 



(38) 



In four dimensions we have A^i > and A^2 > when we take, for example, C = 0.1, A ~ 0.05 and 0.32 < ro < 0.55. 
But when n > 5 the numerators, iVi and iV2, can not be positive at same time. We can see this by doing an analysis 
of the polynomial N2. Using ([57)1 and ([55]) we obtain 



N2 < -2Ar^+\AC{n - 1) - C{n'^ - 1)) - Crl{n - 2){n - - n)C{n - l)/2 - C{n - 1)) 

= C(n - l)(n - 3)[4Ar^'+i - Cr^{n -l){n~ 2){n ~ Zf\l2 < 0. 



(39) 



So, for n > 5 and A > there is no stable circular orbit. It is interesting to note that the addition of a negative 
cosmological constant stabilize, former not stable, circular orbits for any n > 4. 

In the following we analyze the stability of the circular orbits when both A and a are not null. In this case the 
equation of metric function is written as 



r/'-(n-3)(l-/) + 
which has the general solution 

f{r)± = l- 



a(n-4)(n-3)(l-/) 



r 

25 



(2r/'-(n-5)(l-/))=0, 



1 ± 



4a 



(n- l)(?i-2) 



2A 



C 



(40) 



(41) 



where a = (n — 3)(n — 4)a and C is a integration constant. As before, let /+ be the function / when the sign of 
square root is positive and let /_ be the other case. In the limit |q;| — )■ the function /_ reduce to 



/- 



2Ar2 



C 



(n-l)(n-2) (n- l)(n-2)r"-3' 



(42) 



By considering the Gauss-Bonnet term as a perturbation of the Schwarzschild-dS (Schwarzschild-adS) spacetime with 
A > (A < 0), we conclude that C must be positive. 
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When a is small, the asymptotic behavior of /+ is given by 



r 



2 OA ^2 



2Ar^ C 



5 ^ (n- l)(n-2) ^ (n - l)(n - 2)r"-3 ' ^^^^ 

As we can see, here again, there are two families of solutions corresponding to the sign in front of the square root 
in (dH). We define two branches following the same criterion, the GR branch and the non-GR branch. 



A. The GR branch solution (A / 0) 

The functions ly^ and in this case can be written as (|26l) 

^2 ^ rg[(n - l)(n - 2)(1 - V^) + 8Aa - aCr^"(n - 5)] 

a(n- l)(-2rJfV/3r(n-2) + r3C) ' ^ ^ 

^2 ^ B,r^^' + ^3^0-+^ + i32rg + i^irg + gprp^ 

2<5(n - l)(-2r^'\/A:(n - 2) + r3C')r;^[(n - l)(n - 2) + 4a(2A + Cr^-")] 



where 



= 1 + 7 ^ ^ fsA + ^ ) , (46) 



and 



(n- l)(ri-2) 
S4 = C(8Aa + (n - l)(n - 2)){n^ ~ 

B3 = C[-S2ay/l%{n - l)(n - 2) + 2a(n2 _ 8?i + 39)(8aA + (n - l)(n - 2))], 
^2 = 8[[(- l)(n - l)(n - 2) + 8aA][(n - l)(n - 2) + 8aA]] , 
Bi = 2C^a{n - 1)[-V^(5 -n) + n + 3], 

So = 452c2[(n-5)(n-9)]. (47) 

First, when A < it is possible to find stable circular orbits. For example, if a = ±1, A — —1, C ~ 1 then 7 is 
stable when ro > 1.2 {n > 5). Note that /_ = if and only if 

p{r) = 2(-A)r"+'^ + {n - l)(n - 2)r"+^ + {n-l){n- 2)(5r" - Cr^ = 0. (48) 

Analyzing the variation of signs in the coefficients of p, we can see that, if A < then the number of variations is one 
and there is only one positive value f such that p{f) = 0. And for tq > r we have /(rp) > 0. If A > 0, the existence of 
a root of p(r) = depends on the relation between a, A and C. In the cases a = ±1, A = —1, C = 1, the polynomial 
p{r) = 2r"'+'* + (n — l)(n — 2)(r^ — l)r" — has only one positive root. Hence, we have a horizon (at rn < 1) with 
the zone of stability outside this horizon. 

Figure [3] shows the behavior of the functions i'^ and for A — —1, C = 1, a = 1 and 1 < t-q < 3, and n — 5(top 
curve), 6, 7, 8, 9, 10 (bottom curve). We see a small region of instability for large n. Similar behavior is found for the 
same parameters, but a — —1. 

The cases where A > was numerically analyzed. For example, we have stable circular orbits when n = 5, 
5 = C = A = 1 and 0.1 < ro < 0.5; n = 6, a = C = A = 1 and 0.2 < ro < 0.55; n = 7, 5 = 2, C = A = 1 
and 0.3 < tq < 0.55; n = 8, a = 9, C = A = 1 and 0.45 < ro < 0.55, and n = 9, 5 = 100, C = 1, A = 0.05 and 
0.62 < rp < 0.83. In all these examples we have /-(ro) > into the fixed intervals. 

The other cases with A > (n > 10 and a > or n > 5 and a < 0) we calculate the functions ly^ and into a 
grid of points (ro, A, a, n) and at each point of the grid these function are not positive at same time. The grid is built 
with 0.1 < ro < 2 (step 0.1), < A < 2 (step 0.5), 5 < n < 20 and 0.01 < 5 < 100 (step 0.01 into [0.01, 0.1], step 0.1 
into [0.1, 1], step 1 into [1, 100]) and C = 1. In the case a < we work with —100 < a < —0.01, with the same steps 
for the corresponding intervals. 

In these cases (A > 0), each value of n has a different set of parameters (ro, a, C) where both and are positive, 
and this differences turns quit impossible to repeat the graphical analysis as in the A < cases. 
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FIG. 4: Non-GR branch, A 7^ case. The frequencies and for A = 1, C = 1, 3 = 1, 1 < n) < 3, and n = 5, (top curve) 
6, 7, 8, 9, 10 (bottom curve). In this case we have stable orbits of arbitrary size. 



B. The non-GR branch solution (A / 0) 

Again, for completeness we describe the case where /+ is the sohition of ((40|) . In this case, the functions and k} 
(PS)) are given as replacing —^f^ by 

The condition to existence of / (/3* > 0) gives us \{n - l)(n - 2) + 8Aa + 45Cr^-"] > 0. For a < the 
numerator of , = rg [(n — l)(n — 2)(1 + y/]3l) + 8Aa — aCrl~"'{n — 5)], is always positive but the denominator, 
Di, = a{n — l)(+2rg -^^(n — 2) + TqC), is negative, so i'^ < 0, and in this case no circular orbit is stable. 

For a > we first analyze the case A > 0. In this case all parameters are positive and the denominators of z/^ and 

are always positive. The same happens with /+, B4, B3, B2 and the expression [(n — l)(n — 2)(H- ^/JJ^) + 8 Aa] 
(Note that Bi is the coefficient Bi after replacing —\f^ by i = 2, 3, 4). These positive coefficients assure that, 
for To sufficiently large, the numerator of and are positive. Hence, we conclude that there are stable circular 
orbits for n > 5, a > and A > 0. 

If a > and A < 0, it is necessary to set conditions to assure > 0. It happens when [n — l)(n — 2) + SAa > 0, 
what is true when 2|A|a < 3. In this case, from the expressions for and k^, we can see that the same argument 
used in the case A > is valid. Hence, we conclude that there are stable circular orbits for n > 5, a > and A < 0. 

Figure ^ shows the functions and for A = 1, C = 1, a = 1, 1 < rg < 3, and n ~ 5, (top curve) 6, 7, 8, 9, 10 
(bottom curve). In this case we have stable orbits of arbitrary size. 

We find a similar behavior for the equivalent case with A = — 1. 
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Newtonian case 


n = 4 


n = 5 


n > 6 


A = 


SSC 


NSC 


NSC 


A < 


SSC 


SSC 


SSC 


A > 


SSC 


NSC 


NSC 



TABLE I: Stability of circular orbits in n-dimensional Newtonian gravity with or 
without cosmological constant. The existence of some stable curves is denoted by 
SSC, and their not existence by NSC. 



GR case 


n = 4 


n > 5 


A = 


SSC 


NSC 


A < 


SSC 


SSC 


A > 


SSC 


NSC 



TABLE IL Stability of circular orbits in n-dimensional Einstein gravity with or with- 
out cosmological constant. The existence of some stable curves is denoted by SSC, 
and their not existence by NSC. 



VI. SUMMARY AND FINAL REMARKS 



An observationally viable theory of gravitation must allow for the existence of stable limited orbits. It is known 
that, in the Newtonian gravity as well as in classical general relativity, these orbits exist only if the dimension of the 
usual space is three (n = 4). We analyze the linear stability of circular geodesies in a n-dimensional (n > 4) spacetime 
with not null cosmological constant A and we find that there exist cases in which these orbits are stable. 

We also analyse the linear stability of circular geodesies for n-dimensional spherically symmetric static solutions 
of Einstein-Gauss-Bonnet (EGB) gravity. Here we also find that there exist cases with stable circular orbits. The 
differential equation that gives us the metric function has two families of solutions corresponding to the sign in front 
of the square root in ([29]) . Following the definition given at [13], the family with minus (plus) sign is called GR 
branch (non-GR branch) solution. The EGB part of our work is divided in several cases, depending on whether the 
cosmological constant A is zero or not, and on the kind of solution used (GR or non-GR branch solution). 

For the non-GR solution, the sign of the coupling parameter a is very important. When a < 0, no matter the value 
of A or n, it is not possible to find stable circular orbits. Otherwise, when a > 0, for all n > 5, it is possible to find 
stable circular orbits, no matter whether A is null, positive or negative. 

The results obtained for the GR solution are similar. When A is positive or null, and a is negative, the circular 
orbits are not stable. But we show values ofQ;>0, 5<n<9, A>0 and tq such that the existence of stable circular 
orbits is possible. The great difference of this branch solution occurs in the case A < 0. In this case, for example, 
taking A = — l,Q; = ±l,n>5 and rg > 1.2 we have stability. The results above summarized in tables pl lIII[) . 

A similar analysis as the one pres ented in this article is under study for higher dimensions rotating black holes 
like the one presented in Refs. [22| | (five-dimensional GR with magnetic field), [24| (n-dimensional GR) and psj 
(five-dimensional EGB theory). Another possible generalizations and applications of the present work concerns with 
the study of gravitational radiation in higher dimensions in the case where circular orbits do exist. For gravitational 
radiation in n-dimensional spacetimes, see for instance 25| . Also the ultrarelativistic plunge into a multidimensional 
black hole was study in [26j . A generalization of this calculation to circular orbits may be relevant to improve our 
understanding of mini-black hole production in TeV-gravity scenarios. 



EGB 
(a 7^0) 


A = 


A > 


A < 


a > 


a < 


a > 


a < 


a > 


Q < 


non-GR branch 


SSC 


NSC 


SSC 


NSC 


SSC 


NSC 


GR branch 


SSC 


NSC 


SSC 


NSC 


SSC 


SSC 



TABLE III: Stability of circular orbits in Einstein-Gauss-Bonnet theory with or without cosmological constant and different 
signs for the EGB coupling constant. The existence of some stable curves is denoted by SSC, and their not existence by NSC. 
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